The Exact Solution of 1-D SU(n) Hubbard Model 
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In this paper we derive out the exact solution of the SU (n) Hubbard model through the coordinate 
and the algebraic Bethe ansatz methods. The energy spectrum and the Bethe ansatz equations are 
obtained. Furthemore, we analyse the ground state and give out the exact analytic solution of the 
model. 
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I. INTRODUCTION 



The 1-D Hubbard model is one of the significant mod- 
els in the study of the strongly correlated electron sys- 
tems which might reveal an enlightening role in under- 
standing the mysteries of the high-T c superconductivity. 
It also favours a lot of properties of integrable models in 
non-pertubative quantum field theory and mathematical 
physics. It's exact solution was first given by Lieb and 
Wu ji}. Based on the Bethe ansatz equations (Lieb-Wu's 
equations), the 1-D Hubbard model was extensively dis- 
cussed in Refs. Although there are lots of works 
on Hubbard model, the integrability was finished until 
1986 by Shastry @, Olmedilla and Wadati gfj in both 
boson and fermion grade versions. Moreover, the eigen- 
value of the transfer matrix related to the Hubbard model 
was suggested in Rcf. [|l5) and proved through different 
method in Ref. 0,^8| ■ Besides, the integrability and 
the exact solution of the 1-D Hubbard model with open 
boundary condition have been investigated by several au- 
thors, iiy^. 

Recently, based on the Lie algebra knowledge, Maas- 
sarani and Mathieu constructed the Hamiltonian of the 
SU(n) XX model and showed its integrability p^ |. Con- 
sidering two coupled SU (n) XX models, Maassarani suc- 
ceeded in generalizing Shastry's method to construct the 
SU(n) Hubbard model p3j . Furthermore, he found the 
related i?-matrix which ensures the integrability of the 
one-dimensional SU(n) Hubbard model pjfl.flt was also 
proved by Martins for n = 3,4. |p5| , and by Yue and 
Sasaki f26|| for general n in terms of Lax-pair formalism.) 
The exact solution of the SU (3) Hubbard model was also 
given in Ref. J27|. However, the eigenvalue and the eigen- 
vectors of the SU(n) Hubbard model have not been given 
yet. 

In this paper, we apply the same method as done in 
Ref. f27| to find the exact solution of the 1-D SU(n) 
Hubbard model. After recalling some basic notations, 



we write down the wave functions and derive out the 2- 
particlc scattering matrix which governs the amplitude 
of the wave functions in section 2. In section 3 with the 
help of the Yang-Baxter relation we employ the algebraic 
Bethe ansatz method to discuss the amplitude of wave 
functions and the Bethe ansatz equations. The exact 
solution of the SU(n) Hubbard model was then given out. 
Under the thermodynamical limit, the explicitly analytic 
form of the ground state energy is given in section 4. 
In section 5, we make some conclusions and list some 
questions to be considered. 



II. THE COORDINATE BETHE ANSATZ 



The Hamiltonian of the SU(n) Hubbard model is 
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where U is the coulumb coupling constant, and E a 
a, t) is a matrix with zeros everywhere except for an one 
at the intersection of row a and column (5: 
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The subscripts a and k stand for two different E opera- 
tors at site k(k = 1, • • • , L). The n x n diagonal matrix 
C is defined by C = T ^ E aa - E nn . We have also 
assumed the periodic boundary condition E^_ L = E^P . 

It was shown that the Hamiltonian (0) has a (su(n — 
l)©u(l))(7©(su(n— l)©u(l)) r symmetry. The generators 
are 



J. 



a/3 



E 

k=l 



E 



a/3 

CL.k 



(3) 



1 



and 



with 



L 

K a = y^C a ,fc, a, ft 

k=l 



■,n-l, 



a, t. 



(4) 



It is worthy to point out that the system enjoys the 5*0(4) 
symmetry when n — 2 |28f| . But the generator are differ- 
ent from Eq. (|3j) - 

Before proceeding the coordinate Bethe ansatz ap- 
proach, we first begin by introduce some notations as in 
Rcf . [^7| . In SU (n) Hubbard model, there are two types 
of particles named a and t, and each particle can occupy 
(n— 1) possible states. The same type of particles cannot 
appear in one site, but two different types of particles can 
occupy one same site. We denote |n), the vacuum state 



of j-th site, |1) 



l)j the particle states of j-th 



site respectively. Under the appropriate basis 
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we can prove that E" n and E™ a act as a creating and a 
annihilating operators of \a)j state particle respectively. 

Define the particle number operators of the particles 
labeled with (s, a) as: 

L L 



fc=l 



k=l 



(a = 1, 2, • • • , n — 1, s = cr, r) 

we can easily prove that these 2(n — 1) particle number 
operators commutate with the Hamiltonian (m: 



[N?,H}=0, 



(6) 



and commute with each other. So such particle numbers 
are conserved. 

Our main purpose is to give out the eigenvalues and 
the eigenstates of the Hamiltonian ([!]). In this section, 
we first give out the two particle sacttering matrix of the 
system. 

In the coordinate Bethe ansatz method, the eigenstates 
of the Hamiltonian can be assumed as 



(7) 



here |0) is the vacuum state of the total chain and |0) = 
if- 



n.,=i ®\ n )j and 
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xO{x Ql <x Q . 2 <-- - < x QjVo )exp(i> k Pj x Qj ), (8) 



9{xi < x 2 < ■ ■ ■ < Xn ) 




< x No , 
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Here oti — 1, 2, • • • , n — l(i = 1, 2, • • • , No) stands for the 
different particle states, the position of the particle, 
<7i = cr, r the type of ith particle. The summation P and 
Q are taken over all permutations of Nq momenta kj and 
A^o coordinates Xj respectively. The symbols ep and cq 
are the parities of two kinds of permutations. Substitut- 
ing the wave function into the Schrodinger equation 



H\i/j No ) = E\ip No 



(10) 



we can get 



A-' a ^ a i '-(. . . l.. u. - 

^(smfc^sinfc,)^:::;^;^;::^. • •./,,./,•,. • • ■). (n) 

with S"iat (sin ki, sin fcj) being the two-particle scattering 
matrix: 

sin ki — sin kj + i^Poloj 



S^(sink i , S mk j ) = 



sin ki — sin kj + 17 



(12) 



where 7 = ^U^- , Pa-aj is the direct product of two kinds 
of permutation operators which permutes the particle 
styles and particle states simultaneously. The energy of 
the Hamiltonian on this wave function is 



E = 2 cos h + -{L - 2N 



(13) 



It is worthy to point out that the ansatz (Q) is not 
complete ]29| ], which can be observed from a simple ex- 
ample. The state |0) and the states without any \n)j are 
all the eigenstates of the Hamiltonian ([!]) with the same 
eigenenergy E — = \UL, so for this eigenvalue, 
the degeneracy of the eigenstate is {n — l) 2L + 1. 

For the convenience, we denote by Q the amplitude 
A...laj',a'i',--- (• • • , kj, ki, ■ ■ ■). When the j-th particle moves 
across the else particle, it gets an S- matrix Sji(qj,qi). 
Here qj = sin kj . The periodic boundary condition leads 
the following constrains on the amplitude Q 



33+2 ' ' ' 



&3 No Sj 1 Sj2 ■ ■ ■ Sjj- 1 Cj 



e lk J L (j. 
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Eq.([l4[) is similar with the Yang's eigenvalue problem 
fiOf . Its solution will give out the Bethe ansatz equations 
and the amplitude. 
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III. THE EIGENVALUE AND THE 
EIGENSTATES 



In the above section, we have obtained the scattering 
matrix. It can be proved that it satisfies the Yang-Baxter 
equation 30 



Sji(qj,qi)S jk (qj,qk)Sik{qi,qk) = 

Sik(qi,qk)Sjk{qj,qk)Sji(qj,qi). (15) 

Define the monodromy matrix as 

T T (q) = S Tl (q,q 1 )S T2 (q,q 2 ) ■ ■ ■ S TNo (q,q Na ), (16) 

here r is the auxiliary space. We can easily prove that 
the monodromy matrix satisfies the Yang-Baxter relation 

S rr ,{q,q')T T {q)T r ,{q') = T T ,(q')T T (q)S TT ,(q,q'). (17) 

So the eigenvalue problem ( |l4|) is equivalent to finding 
the solution of following problem: 



*r T T T (tfrX = XjC = e{qj)C. 



(18) 



Here e(qj) and £ are the eigenvalue and the eigenstate 
respectively of the transfer matrix Xj which is the matrix 
trace of the monodromy matrix. 



A. The fundamental commutation rules 



From the point of view of a vertex model, we can 
interpret the matrix S T j as the vertex operator, the 
matrix S TT i as the i?-matrix. So the i?-matrix is an 
(2(ra - l)) 2 x (2(n - l)) 2 matrix 



Rji(qj,qi) = a 2 {q^qi) + a 3 (qj,qi)Pji 



(19) 



where 



®2{qj,qi) = — , a 3 (qj,qi) = , — , (20) 



qj -qi+il 



qj -qi+i-y 



and Pji is a (2(n — l)) 2 x (2(n — l)) 2 dimensional permu- 
tation operator. 

The vertex operator is 

£j(<l-<lj) = S T j(q,qj) = 

2(n-l) 

a 2 + a 3 E aP ®E^ a , (21) 

a, 13=1 

where E a/3 is a 2(n— 1) x 2(n — 1) matrix with ev- 
erywhere except for an 1 at the intersection of raw a 
and column (3 and E^ 13 is an operator acting on the j-th 
space. 



The local reference state can be defined as: 

(i) / -i \ (i) 



|0)« = 



(i) 







w 



Tj (n-1) 







-l)xl V ^ / i 2(n-l)xl 



where Sj and Tj are the j-th particle style space and 
particle state space respectively. 

The action of L-operator on this state has the following 
property: 

/ a x a z Ef a,Ef ■ ■ ■ a.Ef^ 1 \ 



(1) 



a 2 






a 2 






a 2 



J 

(22) 



here ax{q i q j ) = a 2 (q,qj) + a 3 (q,qj) = 1. 

The global reference state lO)^ 1 -* is then defined by the 
tensor product of local reference states, i.e. |0)^^ — 
Y\j=i ^lO}^ 1 '- ( Here we use the superscript (1) to distin- 
guish lO)^ 1 ) with |0) appeared in the above section). The 
property of the L-operator suggests the representation of 
the monodromy matrix to take the following form 



Tr(q) 



B{q) B(q) 
C(q) A(q) 



(23) 



2(n-l)x2(n-l) 



where B(g) and C(g) are vectors with dimensions 1 x 
(2n — 3) and (2n — 3) x 1 respectively. The operator A(q) 
is a (2n — 3) x (2n — 3) matrix and we denote its elements 
by A a b(q). The remaining operator B(q) is a scalar. 

In the framework of the above partition the eigenvalue 
problem for the transfer matrix (n~9) becomes 



2n-3 



[B(q)+'£A aa (q)}\$) = A(q)\$). 



(24) 



where A(g) and |<E>) correspond to the eigenvalue and the 
eigenstate respectively. From Eq. (P2[), we know the 
action of the monodromy matrix on the reference state 
|0>«: 



B(g)|0)« = |0)( 1 ) ) C(g)|0}W=0, 



A ab (q)\0)W =S ab l[a 2 (q,q j )\0)W, 

3 = 1 



(25) 



where a,b = 1, 2, ■ • ■ , 2n — 3. The operator B(g) plays the 
role of creation operators over the reference state |0)^^. 



3 



To make further progress we have to recast the Yang- 
Baxter algebra in the form of commutation relations for 
the elements of the monodromy matrix. We shall start 
our discussion by the commutation rule between the op- 
erators B(g) and B(p): 

B(g) ® B(p) = [B(p) <g> B(g)] • f « (q,p). (26) 

where rW(g,p) is an auxiliary (2n — 3) 2 x (2n — 3) 2 matrix 
given by 



^(M 13 , ■ ■ • ,p$) - B(pW) ® Sjvx-x^, • • • ,P$), 

(31) 

here we have formally identified <&o with the unity vector. 
This vector has the following symmetry: 

(n (1) ••• n (1) n (1) ■■■ n {1 ^ - 
9 Ni(Pi >' >Pj-l>Pj >• • • jPjvJ — 

^/VilPi r iPj >Pj-ii iPtfi) r j-ij\Pj-vPj )> 

where the subscripts in ^j_i,-(Pj_i!Pj ) emphasize the 
positions in the iVi-particle space V\ <g> ■ • • ® Vj-\ ®Vj® 
■ ■ ■ ® Vn-l on which this matrix acts non-trivially. Here 
we have already assumed that the (N\ — l)-particle state 
was already symmetrized. 

Applying the diagonal elements of monodromy matrix 
on this eigenstate, we have 

Ni , (1) x 

s( e )i^ 1 ({ P «}))=n^v4i^ 1 (te (i) })) 

i=i ot 2 {Pi ,q) 

N l I (1) \ Hi , (1) 

2^1 , (1) v 11 / (1) (IK 

X|^;_ 1 ( g ,pf ) ;{p| 1) }), (33) 



.(i) ... Jih 



(g,p) = a 3 (q,p) + a 2 {q,p)Pi 



(i) 



(27) 



where P^" 1 is an (2n — 3) 2 x (2n — 3) 2 dimensional per- 
mutation operator. 

The commutation relations between the diagonal and 
creation operator B(g) are 



A(g) ® B(p) = ^4t B (P) ® Ml)} ■ f (1) (Q,P) 



a 2 {q,p) 

«3(g,P) 

OL2{q,p) 



B(g)® A(p), 



(28) 



B(g)B(p) 



gife g) 
a2(p, g) 



B(p)B(g) 



"3(P, g) 

a2(p, g) 



B(g)B(p). 



(29) 



Now we have set up the basic tools to construct the 
eigenstates of the eigenvalue problem (|24|). In the next 
section we will show how this problem can be solved with 
the help of the commutations rules ©, @, ©-(If). 



B. The nested Bethe ansatz 



The purpose of this subsection is to solve the eigenvalue 
problem of the transfer matrix p^). 

First, we construct the eigenstates of the transfer ma- 
trix. The eigenstates of the transfer matrix are in prin- 
ciple built up in terms of a linear combination of the 
products of many creation operators acting on the refer- 
ence state, which are characterized by a set of rapidities 
parameterizing the creation operators. We define the ar- 
bitrary iVi -particle eigenstate as: 

I*^({P« (1) })> = *Ni(p { l\p { 2 l \ ■ ■•,?$)• ^|0} (1) (30) 

where the mathematical structure of the eigenvector 

*JVi (Pi> Pa > ' ' ' >Pni) wm be described in terms of the 
creation operators. We denote the components of vec- 
tor T by jF ai "' ajv i which will be determined later on, 
where the index <ij runs over (2n — 3) possible values 
a, = 1, 2, • • • , 2n — 3. 

We construct the eigenvectors ffrjv^Pi 1 ^, ■ • ■ ,PnI) as 



2n— 3 



N u 



E A aa (q)\^ Nl ({p^})) = ]Ja 2 (q, qj ) 



a=l 



l=1 a 2 {q,p l ) 

Ni , (1), TVo iVi , (1) (1), 

2.— (IK H«a(Pi »©) 11 - (i) (ik 

i=i «2(g,P- ; ) j=i fe=i,fc# i "2(Pi .Pfc ) 



xAW(pW ! {pi 1) })l< ) -i(?,pf ) ,{pF 3 })) ) 
where 

B(g)®$ JVl _ 1 (^ 1) ,...,P l (1) ,...,p«) 



(34) 



(i) 



and 



k=l 



(35) 



(36) 



The symbol p\ means that the rapidity p^ is absent 

from the set {pi , ■ ■ ■ ,Pn]}- 

The terms proportional to the A^-particle eigenstate 
|$jVi ({p! 1 })) are denominated the wanted terms because 
they contribute directly to the eigenvalue. From Eqs. 
(p3|)-(|36|), we can directly get the eigenvalue of the N±- 
particle state 



4 



N l I (1) \ 

A( ?I {^})=n l(Pi ,q) 



i a 2 {p { p,q) 

N JVi , (1), 

n^fe^n^HijT^ferf'w- (37) 

<=i «2(g>P; ) 



3=1 



The remaining ones are called unwanted terms and can 
be eliminated by imposing further restrictions which are 
known as the Bethe ansatz equations: 



n 



gib) .gj) 



n 



c (i) (i)^ 



i= 1,2, •••,Aq. (38) 

In fact, the undeterminated eigenvalue A^(q, {p\ ) }) 
must satisfy the following auxiliary problem 

(39) 

where the inhomogeneous transfer matrix T^\q, {p\ [ }) 
is 

This auxiliary eigenvalue problem has the same structure 
as the eigenvalue problem of the transfer matrix, so we 
can write the solution of this auxiliary eigenvalue prob- 
lem from the foregoing procedure directly. 

The eigenvalue of this auxiliary eigenvalue problem is 

A«(,,^ i) },u 2) })=n ai(pf),9) 



=1 a 2 {pf\q) 



j=i i=i o: 2 {q,Pi ) 

(2) 

and the parameters {p\ } satisfy the following Bethe 
ansatz equations: 

ft ft Mgg.pg) 

xA^( P f,{p«},{p[ 2 )}), 

i = l J ".,JV 2 . (42) 

Here A^ 2 ^(g, {p; }, {pi }) is determined by the second 
auxiliary eigenvalue problem. 

Comparing with above procedure, we know that this 
second auxiliary eigenvalue problem is equivalent to find- 
ing the eigenvalue of the transfer matrix of the SU(n — 1) 
Hubbard model. Applying the foregoing procedure re- 
peatedly until the dimension of the lower level space be- 
comes equal to one, at this point, the final auxiliary r- 
matrix is 



f (2n - 3) (q,p) = l, 



(43) 



and the nested Bethe ansatz closed. For n = 2 there is 
no nested necessary. 

Now we have obtained the eigenvalue of the transfer 
matrix: 



iVi ( (1) 

A( q ,{p ( i 1) }) = U" " 



(i) u _ TT a i(Pl >g) 
i at 2 {Pi ,q) 



+ n *> n ^Hi)T A(i) (^ ( 44 ) 

i=i «2(g,P- ') 



3=1 



and 



N„ 



A (m) (qApr}Apr +1} })= n 



*i(p { k m+1 \q) 
ijl / (m+1) \ 

fe =l a2(pL ',g) 



(m)x TT «l(g,Pfc 



(m+l)s. 



i=i fc=i "2(g,^ ') 

KA^'fe^ 1 '}), (45) 
where m = 1, 2, • • • , 2n — 3. We have defined N 2n - 2 = 0, 
and A( 2 "-%, {pf n " 2) }) - 1. 

We also have the following (2n — 3) sets of Bethe ansatz 
equations: 



N„ 

n 



a HPi 'P.7 ) TT «l(P; _>Pj ) 



n 



lj. / (m+1) (m)s 11 /Jm+2) 

=1 «2(Pi ,P» ) j =1 a 2 (Pi -Pi J 



n 



»+l / (m+1) (m+1) 

; ,-l X-Li ^(Pj ■ 



J™+1) j™+ih 



(46) 



where to = 1, 2, • • • , 2n — 3 and the periodic boundary 
condition (|l4|) implies: 



JVi / (1) \ 

ik s L _ TT "HP; .gj) 

1=1 MpI ,qj) 



(47) 



so we get 2(n — 1) sets of Bethe ansatz equations. 

After Substituting the expressions of ot\, a 2 into above 
expressions, the 2(n — 1) sets of Bethe ansatz equations 
become to be 



and 



n 



ft (pf-p, (i) -»i) 

l\^-P^ + qy 

, -im+1) ~(m+l) . ■ \ 

(Pi - Pfc + h) _ 

--(m+1) -(m+1) 



(48) 



Wm+l f~(m+l) _ -(m+1) + 
k=l,kj&i (Pi Pk 1 I) 



M(pi m+1) -^ m) -^) 

2V„ 



n 



± 2 (f>t +1) -P { r +2) + q) 



(m = 1,2, • • • ,2n - 3). 



(49) 
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Here p^ = sin kj and we have used the shifted param- 
eters pj = p\ k ^ + ik-k to bring our equations to more 
symmetric forms. 

In principle, we can consider the same system with 
twisted boundary condition, similar as done in SU (2) 
Hubbard model. The Bethe ansatz equations and the 
energy will pick up an appropriate factor which related 
to a twisted angles. On the other hand, we can also 
study this method based on the fermionic base, which 
can be obtained by using the Jordan- Wigner transfor- 
mation from bosonic one, but the boundary condition 
will be changed. Generally, the bosonic SU(n) model 
with periodic boundary will correspond to a fermionic 
one with twisted boundary and viceversa, the fermionic 
system in the presence of periodic boundary correspond 
to a twisted bosonic version. This is a direct generaliza- 
tion of discussions give in Ref. O] 



IV. THE GROUND STATE 



In the above section, we have obtained the final 2(n— 1) 
sets of Bethe ansatz equations of the SU(n) Hubbard 
model. In this section, we will use them to analyse the 
ground state of the model. 

Under thermodynamic limits, after taking the lo- 
grarithm, the Bethe ansatz equations (^8|)-(^9|) for the 
ground state change into: 



27Tp (0) (fc) = 1-COsfc / 

J-BW 



4 7/0 M(A') 



, 7 2 + 4(A(°) - A') 



rdA', 



(50) 



and 



_BC») 7 2 + (AW - A') 2 
' 47p( m - 1 >(A , ) 
_bc»-d 7 2 + 4(A(") - A') 2 
s<m+1) 4 7P (™+ 1 )(A') 



dA' 



-BCm+D 7 2 + 4(A(") - A') 2 ' 
(m = 0,l,---,2n-3), (51) 



where A^ ) = sin k and i?( m ) are determined by the con- 
ditions 



/ pW(A)dA=^, 



(52) 



and B^-V = 0. The functions p( m )(A) are the dis- 
tribution functions of real parameters kj and p[^ , (i — 
1, 2, • • • , 2n — 3) respectively. 



Eqs. (pOJ)-(plj) have a unique solution which is pos- 



itive for all allowed B^ m \ 



N„ 



,(m = 0,l,---,2n - 3) 



are monotonically function of I?( m ) respectively. Thus 
the ground state is characterized by B^ = tt, i?( m ) = 
oo, (m = 1, 2, • • • , 2n — 3). 

After taking Fourier transforms of the above equations 
we can obtain the result of the distribution functions 



3<°>(fc) = 



1 cos k 
2tt 2tT 



, ^,7 8 '^ M ^^, (53) 



sinh((n — l)7|w|) 



and 



p( m '(A) 
1 

2^ 



sinh((2n - 2 



m ) 2 M 



sinh((n — l)7|o;|) 
(to = 1,2, • ■ • ,2n - 3) 



■J (u)e- tulA duj, 



(54) 



From Eqs.(|52|)-(|54|) we have Nq — L, this means that 
all lattice site is filled by one particles, and N m — 
2 2(n-^L)" ^o , these corresponding to the 2(n— 1) conserved 
particle numbers. 

The ground state energy then is given by 

E = 2L ( p(k) cos kdk - ^-L 

J-OG ^ 

sinh((2n — 3)^|w|)e~^ 



-2L 



wsinh((n — l) 7 |o;|) 



- Jq(oS) J±(u>)du> 
(55) 



V. CONCLUSIONS 



The main purpose of this paper is to investigate the 
eigenvalues and the eigenstates of the Hamiltonian of the 
SU(n) Hubbard model. We have succeeded in finding the 
eigenvalues of the Hamiltonian (0) and obtained 2(n — 
1) sets of Bethe ansatz equations. Based on the Bethe 
ansatz equations, we have found the explicit expression of 
the energy and the distribution functions of the rapidities 
corresponding to the ground state for positive U. 

From the result of present paper, we know that the 
result in Ref. |27]] is not complete. The solution of the 
SU(3) Hubbard model obtained in Ref. |27|] was not the 
ground state solution. Next, from present paper we also 
know that the number of Bethe ansatz equations should 
be equal to the number of the conserved particle num- 
bers. But in Ref. |^7| one set of Bethe ansatz equations 
is missing. 

For SU (n) Hubbard model, another important ques- 
tion is to study the excitation spectrum and the low tem- 
pareture thermodynamics of the model for both positive 



G 



and negative U. The Bethe ansatz equations given in 
present paper will play an key role. As we know, the 
negative U has a distinguished properties from positive 
U case. The solution structure of Bethe ansatz equations 
corresponding to the ground state and excited spectrum 
for attractive case (U < 0) are different from those for 
repulsive case. It is worthy to be studied in future. 

In present paper, we have given energy spectrum of the 
Hamiltonian (|l|). As an integrable model, it is also impor- 
tant to find the eigenvalue of infinit number of conserved 
laws. The useful approach is algebraic Bethe Ansatz 
method. 

In Ref. (24) , the I/-operator and i?-matrix were given, 
therefore, one can define the transfer matrix as done in 
usual Hubbard model (SU(2) case). But how to find the 
eigenvalue of the transfer matrix is unknown. It may be 
solved by using the method proposed in Ref. |l8| . We 
will consider this problem late. 



[21] M. Shiroishi and M. Wadati, J. Phys. Soc. Jpn. 65 (1996) 

1983; J. Phys. Soc. Jpn. 66 (1997) 1. 
[22] Z. Maassarani and P. Mathieu, Nucl. Phys. B 517 (1998) 

395-408. 

[23] Z. Maassarani, Phys. Lett. A 239 (1998) 187-190. 
[24] Z. Maassarani, Mod. Phys. Lett. B 12 (1998) 51-56. 
[25] M. J. Martins, Phys. Lett. A 247 (1998) 218. 
[26] R. Yue and R. Sasaki, P. Phys. Soc. Jpn. 67 (1998) 2967 
[27] B. Y. Hou, D. T. Peng and R. H. Yue, Nucl. Phys. B 575 
(2000) 561 

[28] C. N. Yang and S. Zhang, Mod. Phys. Lett. B 4 (1990) 
759. 

[29] F. H. L. E/31er, V. E. Korepin and K. Schoutens, Nucl. 

Phys. B 384 (1992) 431 
[30] C. N. Yang, Phys. Rev. Lett. 19, (1967) 1372 



E. Lieb and F. Y. Wu, Phys. Rev. Lett. 20 (1967) 1445. 
A. A. Ovchinnikov, Sov. Phys. JETP 30 (1970) 1660. 

G. V. Uimin and S. V. Fomichev, Sov. Phys. JETP 36 

(1973) 1001. 

C. F. Coll, Phys. Rev. B 39 (1974) 2150; 

F. Woynarovich, J. Phys. C 15 (1982) 85, 97, C 16 (1983) 
5293. 

A. KKimpor, A. Schadschneider, J. Zittartz, Z. Phys. 
B78 (1990) 99. 

M. Takahashi, Prog. Theor. Phys. 42 (1969) 4089; 43 
(1970) 860, 1629. 

H. Shiba, Phys. Rev. B 6 (1972) 930. 

F. Woynarovich and K. Penc, Z. Phys. B 85 (1991) 261. 
M. Takahashi, Prog. Theor. Phys. 47 (1972) 69, 52 

(1974) 103. 

T. Kowa, Prog. Theor. Phys. (1990) 655. 

N. Kawakami, T. Usuki, A. Okiji, Phys. Lett. A 137 

(1989) 287. 

H. Frahm, V. E. Korepin, Phys. Rev B 42 (1990) 10573, 
B 43 (1991) 5653. 

J. Carmelo, P. Horsch, A. A. Ovchinnikov, Phys. Rev. B 
45 (1992) 7899. 

B. S. Shastry, Phys. Rev. Lett. 56 (1986) 1529; Phys. 
Rev. Lett. 56 (1986) 2453; J. Stat. Phys. 50 (1988) 57. 
E. Olmedilla and M. Wadati, Phys. Rev. Lett. 60 (1987) 
1595. 

R. H. Yue and T. Deguchi, J. Phys. A: Math gen. 30 

(1997) 849-865. 

M. J. Martins, P. B. Ramos, Nucl. Phys. B 522 [FS] 

(1998) 413-470. 

H. Asakawa and M. Suzuki, J. Phys. A 29 (1995) 225. 
T. Deguchi and R. H. Yue, Exact Solution of 1-D Hubbard 
Model with Open Boundary Conditions and the Confor- 
mal Dimensio ns under Boundary Magnetic Fields (1996) 



OCHAPP-84, cond-mat/9704138 



7 



